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ABSTRACT 

We present a broad-band analysis of deep Suzaku observations of nearby Seyfert 1 
AGN: Fairall 9, MCG-6-30-15, NGC 3516, NGC 3783 and NGC 4051. The use of 
deep observations (exposures > 200 ks) with high S/N ahows the complex spectra of 
these objects to be examined in full, taking into account features such as the soft ex- 
cess, reflection continuum and complex absorption components. After a self-consistent 
modelling of the broad-band data (0.6-100.0 keV, also making use of BAT data from 
Swift) , the subtle curvature which may be introduced as a consequence of warm ab- 
sorbers has a measured affect upon the spectrum at energies > 3keV and the FeK 
region. Forming a model (including absorption) of these AGN allows the true extent 
to which broadened diskline emission is present to be examined and as a result the 
measurement of accretion disc and black hole parameters which are consistent over 
the full 0.6-100.0 keV energy range. 

Fitting relativistic line emission models appear to rule out the presence of max- 
imally spinning black holes in all objects at the 90% confidence level, in particular 
MCG-6-30-15 at > 99.5% confidence. Relativistic FeK line emission is only marginally 
required in NGC 3516 and not required in NGC 4051, over the full energy band- 
pass. Nonetheless, statistically significant broadened 6.4 keV FeKa emission is de- 
tected in Fairall 9, MCG-6-30-15 and NGC 3783 yielding black hole spin estimates 
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and a < —0.04 respectively, when fitted with disc 
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> 1 INTRODUCTION 



1 



Gaining information regarding the spin of supermassive 
black holes (SMBHs) in active galactic nuclei (AGN) is 
essential in order to understand some of the fundamental 
mechanisms for powering radio jets, galaxy evolution and 
AGN merger histories. For example, it is thought that the 
spin of the SMBH may be related to powering relativistic 
jets and it is therefore important to begin to assess the dif- 
ferences between radio-loud and radio-quiet AGN (Wilson & 
Colbert 1995; Moderski et al. 1998; Mangalam et al. 2009) 
and the intrinsic properties of the SMBH and its accretion 
disc. Similarly, the distribution of SMBH spin or the 'av- 
erage' spin could be an essential tool for distinguishing be- 
tween various galaxy evolution models (Hughes & Blandford 
2003; Volonteri et al. 2005; King & Pringle 2007; RezzoUa et 
al. 2008). Prolonged accretion could result in a large fraction 
of BHs having a high degree of spin, however, if the major- 
ity of AGN have more intermediate or low BH spin then 



spin- alignment during mergers could be the mechanism re- 
sponsible for spin-ups or the Blandford-Znajek effect (the 
magnetic extraction of BH rotational energy; Blandford & 
Znajek 1977) could cause the reduction of BH spin in many 
AGN (Berti & Volonteri 2008). 

The properties of the SMBH and the surrounding accre- 
tion disc can be measured through the analysis of the X-ray 
spectrum of AGN since line emission occurring from close to 
the black hole can become relativistically broadened (Fabian 
et al. 1989; Laor 1991). The shape of such a line profile al- 
lows properties such as the emissivity and the inclination 
of the disc to be determined, in addition to the typical in- 
ner radius of emission and the spin of the central black hole 
(Dovciak et al. 2004; Brenneman & Reynolds 2006; Dauser 
et al. 2010). These signatures of emission from close to the 
black hole are most prominent in the neutral FeKa line at 
6.4 keV due to the high abundance and fluorescent yield of 
iron. In order to make the required measurements (partic- 
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ularly BH spin), a robust broad-band spectral model must 
be formed in order to assess the contribution (if any) of rel- 
ativistically broadened line emission. 

Recent publications have gone some way towards mak- 
ing tentative spin measurements in AGN through the analy- 
sis of Seyferts lacking significant amounts of intrinsic absorp- 
tion (Miniutti et al. 2009; SchmoU et al. 2009; Patrick et al. 
2011), or the analysis of deep observations of more complex 
AGN (Brenneman & Reynolds 2006; Brenneman et al. 2011; 
Gallo et al. 2011). As noted in Patrick et al. (2011), the in- 
terpretation as to the origin of the soft excess can influence 
the obtained accretion disc parameters and BH spin. Em- 
ploying a relativistically smeared reflection spectrum from 
the inner regions of the accretion disc (e.g. an ionized disc 
reflection spectrum convolved with a Kerr metric; Ross & 
Fabian 2005; Brenneman & Reynolds 2006) to model both 
the soft excess and broadening in the FeK region yields re- 
sults in which the measured parameters differ significantly 
when they are obtained from either the soft excess or FeK 
region. This is particularly notable in AGN with a strong 
soft excess, in which case the soft excess is the driving force 
behind the fit, often leading to typically high spin and emis- 
sivity values i.e. a > 0.90 and q > 4.5. As the soft excess is 
often smooth or featureless, a very high degree of blurring 
is required. 

Alternatively employing a Comptonization origin of the 
soft excess (Titarchuk 1994) and seperately measuring these 
parameters purely from the observed line profile in the Fe K 
region yields a more moderate emissivity index and a range 
of spin values akin to the results found in objects without an 
observed soft excess i.e. < a < 0.998 and q ~ 2.3 (Patrick 
et al. 2011). It is therefore essential that high quality data 
sets and deep exposures are used in order to distinguish 
between these scenarios and measure the true accretion disc 
parameters independently. 

This paper includes a small sample of AGN, all fea- 
turing high quality, long exposure (> 200 ks) observations 
with Suzaku which are publically available. This is further 
to previous work on the analysis of a sample of relatively 
'bare' AGN (Patrick et al. 2011) applied to this sample of 
more complex AGN i.e. those featuring complex warm ab- 
sorbers. Suzaku is unique in that it offers the opportunity to 
gather both soft and hard X-ray data simultaneously using 
the XIS (Koyama et al. 2007) and HXD detectors (Taka- 
hashi et al. 2007) respectively. This enables full broad-band 
modelling of the 0.6-100.0 keV spectrum, including features 
such as the Compton refiection hump peaking at ~ 30.0 keV 
(George & Fabian 1991) which would otherwise be absent 
in spectra obtained from observations which only cover the 
lower energy range. Only with high energy X-ray data can 
the refiection component be fit simultaneous to the Fe K line, 
allowing its strength and ionization state to be determined 
and subsequently its contribution to the FeK region. Also 
including time averaged BAT data (over a 58 month dura- 
tion containing the Suzaku observations) from the Swift 58 
month all sky survey (Baumgartner et al. 2010) increases 
the energy range up to ~ 100 keV. 

Alternative explanations of apparent 'broad' residuals 
in the Fe K region have also been forward by some analyses 
in the literature, see Turner & Miller (2009) for a review. 
Models consisting of complex warm absorbers and partial 
covering scenarios can often reduce the observed strength of 



any broadened residuals in the Fe K region (e.g. through con- 
tinuum curvature imparted by bound-free transitions) and 
in some cases may account for these entirely. In particular, 
some of the objects in this sample from past observations 
have been modelled in such a way: MCG-6-30-15 (Miller et 
al. 2009), NGC 3516 (Turner et al. 2005) and NGC 3783 
(Reeves et al. 2004). The complex absorption used to model 
these AGN is suggested to reduce the need of additional rel- 
ativistic emission and indeed leads to the question whether 
diskline emission can be seen which is in stark constrast to 
typical refiection dominated models such as those used by 
Miniutti et al. (2007) and Brenneman & Reynolds (2006). 

In this paper, we first fit the spectrum at soft X-ray 
energies and include any absorbing components which may 
be required by the data prior to assessing the remaining 
residuals in the FeK region which may be representative 
of diskline emission. By taking a balanced approach over a 
broad bandpass we aim to determine the effects of warm 
absorbers upon the FeK region in an attempt to discover 
if these broad residuals remain. Only after modelling for 
absorbers do we include any models which may be repre- 
sentative of emission from the inner regions of the accretion 
disc. 

This paper therefore aims to assess the strength or 
indeed even the presence of diskline emission in the FeK 
regions of these objects. Through modelling of absorption 
zones where required and the application of a variety of 
models and interpretations, the degeneracies involved with 
measuring such parameteres can be investigated and hence 
estimates can be placed upon BH spin and accretion disc 
emissivity/ inclination. We perform an analysis of these ob- 
jects by forming models both with and without a partial 
covering component in order to assess the effect of the par- 
tial covering absorber on the iron line parameters. In later 
models only after modelling the absorption do we model 
with a dual inner (blurred) and outer (unblurred) reflector 
from distant material in an attempt to estimate accretion 
disc and BH parameters. 



2 OBSERVATIONS & DATA REDUCTION 
2.1 Observations 

The objects included within this sample are listed in Table 
[T]and are all the Seyfert 1, radio quiet AGN with exposures 
> 200 ks which have been observed with Suzaku with data 
publically available in the Suzaku data archive^. The objects 
in this sample are also all nearby with redshift z < 0.05. 
High energy X-ray data from Swift/BAT is also used in 
addition to that obtained from the HXD detector onboard 
Suzaku (allowing the cross-normalization to vary), therefore 
the total energy range covered is 0.6-100.0 keV. Details of 
the observations included are listed in Table (2] NGC 4051 
has been observed on three occasions with Suzaku, only the 
long November 2008 observation is included here, Lobban et 
al. (2011) conducted an analysis of all three simultaneously. 
The long 2010 Suzaku observation of NGC 1365 is excluded 
from this paper due to its similarities with Seyfert 2 AGN, 



^ http://heasarc.gsfc.nasa.gov/ 
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Table 1. The Suzaku Seyfert sample 



Ubject 


KA (JzUUUj 


Dec (JzUOUj 


Heasniit 


N u \yjci\) (10 cm ) 


Fairall 9 


01 23 45.8 


-58 48 21 


0.0470 


0.0316 


MCG-6-30-15 


13 35 53.8 


-34 17 44 


0.0077 


0.0392 


NGC 3516 


11 06 47.5 


+72 34 07 


0.0088 


0.0345 


NGC 3783 


11 39 01.7 


-37 44 19 


0.0097 


0.0991 


NGC 4051 


12 03 09.6 


+44 31 53 


0.0023 


0.0115 



Table 2. Summary of observations for the objects in the sample. ^ The observed 2— lOkeV flux for XIS, 15-50 keV flux for HXD and 
20-100 kcV flux for BAT, in units lO^'^^ergcm"^ s"'^ from the baseline model. The XIS count rates listed are per XIS. We use BAT data 
from the 58 month BAT catalogue (Baumgartner et al. 2010). 



Object 


Mission 


Instrument 


Date 


Exposure (s) 


Count rate 




Fluxl 


Obs. ID 




SuzaJcu 


XIS 
HXD 


2010/08/26 


^29290 


"1 'yori 1 ri nrio 

i.ToU ± U.UUz 
u.uoy zn. u.uuz 




2.45 

O. Ol 


705063010 




Swift 


BAT 






(4.7 + 0.2) X 10- 


4 


3.34 






SuzaJcu 


XIS 
HXD 


2006/01/09 


143196 
118900 


2.862 ± 0.003 
0.094 ± 0.002 




4.38 
4.70 


700007010 


MCG-6-30-15 


SuzaJcu 


XIS 
HXD 


2006/01/23 


98483 
76800 


2.461 ± 0.003 
0.103 ± 0.003 




3.81 
5.06 


700007020 




SuzaJcu 


XIS 
HXD 


2006/01/27 


96691 
83660 


2.708 ± 0.003 
0.104 + 0.002 




4.16 
4.98 


700007030 




Swift 


BAT 






(6.8 + 0.3) X 10" 


-4 


3.86 




NGC 3516 


SuzaJcu 


XIS 
HXD 


2009/10/28 


251356 
178200 


0.456 + 0.001 
0.059 + 0.001 




1.35 
3.32 


704062010 




Swift 


BAT 






(10.5 + 0.2) X 10- 


-4 


7.76 




NGC 3783 


SuzaJcu 


XIS 
HXD 


2009/07/10 


209503 
162000 


1.942 + 0.002 
0.235 + 0.002 




5.92 
12.04 


704063010 




Swift 


BAT 






(16.6 + 0.2) X 10" 


-4 


11.54 




NGC 4051 


SuzaJcu 


XIS 
HXD 


2008/11/06 


274350 
204500 


1.858 + 0.002 
0.062 + 0.001 




2.45 
3.05 


703023010 




Swift 


BAT 






(4.0 + 0.2) X 10- 


-4 


2.80 





i.e. due to its high degree of obscuration (Maiohno et al. 
2010). 

2.2 Data reduction 

The Suzaku data in this paper were reduced using version 6.8 
of the IfEASOFT data reduction and analysis package. The 
XIS source spectra were extracted using 3.0' circular regions 
centred on the source. Background spectra were also ex- 
tracted using 3.0' circular regions, this time centred on a re- 
gion of the CCD not featuring any of the source or Fe 55 cali- 
bration regions. Only data from the front-illuminated XIS 0, 
2 and 3 detectors were used due to their greater sensitivity 
at FeK energies, however data from the back-illuminated 
XIS 1 remains consistent with the front-illuminated detec- 
tors. Observations since November 2006, however, do not 
include data from XIS 2 as it is no longer operational. 

MCG-6-30-15 was observed on three occasions in Jan- 
uary 2006, this paper makes use of all three observations 
and the time averaged spectrum is used in the main analy- 
sis. Analysis of the individual observations yields FeK line 
profiles which were consistent within errors and as such all 
three of the January 2006 observations have been coadded. 
The events files for each XISO observation (i.e. all three of 



the observations during January 2006) were loaded into xs- 
ELECT and the source and background spectra produced. 
This was also performed for XIS 2 and XIS 3 to produce 
single time-averaged source and background spectra for the 
observations of MCG-6-30-15 listed in Tabled 

The HEASOFT tool xisrmfgen was used to generate 
the XIS redistribution matrix files (rmf) and the ancilliary 
response files were generated using xissimarfgen. The data 
from each of the front-illuminated XIS detectors were then 
co-added using mathpha, ADDRMF and ADDARF in order 
to increase signal to noise. We ignore all XIS data below 
0.6 keV, above 10.0 keV and between 1.7-2.0 keV due to cal- 
ibration uncertainties of the detectors around the Si K edge. 

Spectra from the HXD were extracted from the cleaned 
HXD/PIN events files and subsequently corrected for in- 
strument dead time using the tool hxddtcor. The tuned 
HXD/PIN non X-ray background (NXB) events files were 
used for background subtraction (Fukazawa et al. 2009) and 
generated with 10 x the actual background rate in order to 
reduce photon noise, also with identical good time intervals 
(GTIs) as used in the source events. A simulated cosmic 
X-ray background was also produced for each observation 
using XSPEC v 12.6.0q with a spectral form identical to 
that used in Gruber et al. (1999), this was then added to 
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the corrected NXB to form a single background file for each 
observation. The appropriate response and flat field files 
were used from the Suzaku CALDB suitable for the respec- 
tive epochs of each observation. IfXD data is used over the 
15.0-60.0 keV range along with the appropriate PIN/XIS- 
Fl cross-normalisation according to the epoch and nomi- 
nal pointing position^. Additional hard X-ray data from 
the BAT instrument onboard Swift obtained from the 58 
month BAT catalogue were also included for all objects 
over the range 20.0-100.0 keV with a BAT/XIS-FI cross- 
normalisation allowed to vary (but typically ~ 1). 



3 ANALYSIS & RESULTS 

Spectral analysis and model fitting is performed from within 
XSPEC V 12.6.0q (Arnaud 1996), all models are modified 
by Galactic absorption which is accounted for by the WABS 
multiplicative model (Morrison & McCammon 1983). The 
respective Galactic column densities were obtained using 
the NH ftool for each source giving the weighted average 
A'h value of the LAB Survey of Galactic H I (Kalberla et al. 
2005), using abundances from Anders & Grevesse (1989). 
Data is fit over the full 0.6-100.0 keV range available, ex- 
cluding those regions affected by the uncertainties in the 
XIS calibration mentioned above. The minimisation tech- 
nique is used throughout, all errors are quoted at the 90% 
confidence level (Ax^ ~ 2.71 for one interesting parame- 
ter) and include statistical and not instrumental systematic 
errors. Where the significance of components is quoted in 
terms of Ax^, the component in question has been removed 
from the model and the data refit to ensure that the order 
in which components are added to the model does not play 
a role in the quoted statistical significance. Throughout this 
paper, a positive Ax^ corresponds to a worsening in the fit, 
whereas a negative Ax^ corresponds to an improvement in 
the fit. 

3.1 Spectral Models 

3.1.1 Modelling of the soft excess 

There are a number of possible origins of the soft excess 
observed in AGN X-ray spectra. A simple method of mod- 
elling the soft excess is as a black body to replicate direct 
thermal emission from the accretion disc (Malkan & Sar- 
gent 1982). However, the relatively constant temperature 
of the soft excess is in disagreement with the typical ac- 
cretion disc properties for a SMBH which should scale with 
Mgj^^* - more typical of an intermediate mass black hole e.g. 
kT ~ 100 - 200 eV (Gierlihski & Done 2004). A variation 
upon this concept is that the soft excess originates from the 
inverse-Compton scattering of the EUV photons from the 
accretion disc in a hot plasma above the disc, which would 
allow for a large variation in photon seed temperature, yet 
it would still produce a relatively constant photon temper- 
ature such as the compTT model (Titarchuk 1994). Such a 
model has been successful in modelling the soft excess in a 
sample of PG quasars (see Porquet et al. 2004). 

^ http: / /heasarc . gsfc . nasa.gov /docs/ suzaku / analysis / watchout .html 



Departing from an origin of the soft excess due to emis- 
sion from the accretion disk, an alternative explanation is 
an atomic origin. In this scenario soft X-ray emission lines 
are responsible for the soft excess, however since these are 
discrete features and the observed soft excess is typically 
smooth, the refiection continuum must be relativistically 
blurred and hence originate from the inner regions of the 
accretion disc (Ross, Fabian & Ballantyne 2002). This in- 
terpretation of the soft excess has also proved successful in 
some AGN such as Ark 120, Fairall 9 and MCG-6-30-15 
(Ballantyne, Vaughan & Fabian 2003; Crummy et al. 2006; 
SchmoU et al. 2009; Nardini et al. 2010). 

An alternative origin for the soft excess as suggested 
by Gierlihski & Done (2004) is relativistically smeared ab- 
sorption whereby the discrete atomic features are masked by 
high velocities. This would result in a 'hole' in the spectrum 
creating an observed soft excess and hardening at higher 
energies. 

3.1.2 Compton reflection 

Following the availability of hard X-ray data from Suzaku's 
HXD and Swift's BAT detectors features such as the Comp- 
ton refiection hump at ~ 30keV can be measured in the 
spectra of these AGN. Using an ionized reflection model such 
as REFLIONX (Ross & Fabian 2005) allows properties such 
as the ionization state, the iron abundance (Zpc) of the disc 
and strength of the refiection component to be measured. 

Narrow fiuorescent emission lines in AGN spectra result 
from refiection off distant material and as such an unblurred 
refiection component should be included as part of a base- 
line model in AGN spectra. Previous studies of the iron line 
regions in Seyfert AGN by Bianchi et al. (2004), Nandra et 
al. (2007) and Patrick et al. (2011) also suggest that ionized 
species of iron are relatively common at energies of 6.7 keV 
and 6.97 keV for FeXXV and FeXXVI respectively. Neu- 
tral narrow Fe Ka and the accompanying Fe K/3 emission at 
6.4 keV and 7.056 keV respectively are ubiquitous in AGN 
spectra (Nandra et al. 2007). 

To account for the distant refiection component in these 
AGN, we use the pexrav (Magdziarz & Zdziarski 1995) 
neutral refiection model initially which is then replaced in 
later models with an unblurred reflionx ionized reflection 
model with the input photon index F tied to that of the in- 
trinsic powerlaw. The reflionx model includes soft narrow 
emission lines in addition to narrow FeKa. However, FeK/3 
emission with fiux Fkp ~ 0.13 x Fkc is not included self- 
consistently in reflionx and as such is modelled using a 
narrow Gaussian of fixed width a = 0.01 keV and with fiux 
fixed at the value obtained during the initial parametrisa- 
tion. Additional narrow ionized emission lines e.g. at 6.7 keV 
and 6.97 keV may also be required and are modelled using 
a narrow Gaussian with width fixed as above. 

3.1.3 Warm absorption 

The presence of an X-ray absorber primarily affects the X- 
ray spectrum of soft energies, however, with higher column 
densities the warm absorber may add subtle spectral cur- 
vature to the spectrum above 2.5 keV. In particular, highly 
ionized absorption zones can affect the measured line pa- 
rameters and strength of the observed broad red-wing in 
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the FeK region (Reeves et al. 2004; Turner at al. 2005). 
In previous observations of these objects (bar Fairall 9) the 
presence of complex absorption zones has been well docu- 
mented (Kaspi et al. 2000; Yaqoob et al. 2005; McKernan et 
al. 2007; Markowitz et al. 2008; Miller et al. 2008; Lobban 
et al. 2011). 

To model the soft X-ray warm absorber components in 
this paper, we use the xstar generated (Kallman et al. 2004) 
grid illuminated by a photon index of F = 2, abundances 
fixed at solar values (except Ni, which was set to zero) and 
with a turbulent velocity of 200kms~^. This grid is well 
suited to accounting for typical absorption zones due to its 
wide range in column density (5 x lO^^cm"^ < TVh < 5 x 
lO'^'^cm"^) and ionization parameter (0 < log^ < 5). During 
the fitting process absorption zones are added as required, 
in some objects more than one zone may be statistically 
required. 



3.1.4 Highly ionized absorption 

In the event that absorption features in the FeK region 
are found (such as ls-2p resonance lines from FeXXV or 
Fe XXVI) , these will be accounted for using a model repre- 
sentative of a highly ionized absorption zone. We use an XS- 
TAR generated grid with a turbulent velocity of 1000 km s~^ 
(to better model the observed equivalent widths of the lines), 
an input continuum of F — 2.0 and range in column den- 
sity and ionization parameter of 1 x lO^'^cm"^ < A'^e < 
1 x 10^*cm~^ and < log^ < 6 respectively. 



3.1.5 Partial covering 

Some models used in the analysis of these AGN use partial 
covering geometries whereby a fraction of the observed X- 
rays are absorbed by a surrounding gas in the line-of-sight 
(in addition to typical fully-covering absorbers) as described 
in Section 3.1.3, while some fraction of the continuum 'leaks' 
through and is unaffected by the partially covering material. 
The partial coverer used here takes the form of (powerlaw 
-)- WA*powerlaw) with the photon index of both the pow- 
erlaw components tied and normalisation free to vary. The 
parameters of the warm absorber (ionization and column 
density) are also allowed to vary. The column density of 
the partially covering medium can have a significant affect 
upon the spectrum. For example, high column partial cover- 
ers (A'h ~ lO'^* cm~'^) predominantly affect the hard X-ray 
energies and are used here to supplement the distant reflec- 
tion component in cases where a hard excess remains. Lower 
column density partial coverers may greatly affect spectral 
curvature at lower X-ray energies and in some cases can en- 
tirely remove any 'broad' residuals in the FeK region (Miller 
et al. 2009). 



3.2 Initial parametrization and baseline model 

An initial parametrization of the spectra is made in order 
to assess properties such as the strength of the reflection 
continuum and the strength of any emission or absorption 
features which may be present in the Fe K region. It should 
be noted that all objects in this sample (with the exclusion 



of Fairall 9) feature complex intrinsic absorption compo- 
nents. Once these basic components have been established, 
resulting in a reasonable fit to the data, the Fe K region can 
then be assessed as to the presence of any additional fea- 
tures. For example, narrow emission lines such as FeXXV 
and FeXXVI, or indeed the need for a highly ionized ab- 
sorption zone (see Reeves et al. 2004 for Fe XXV absorption 
in an XMM-Newton observation of NGC 3783) may be re- 
quired to model the Fe K region. 

The strength of the reflection component is ini- 
tially measured using the pexrav neutral reflection model 
(Magdziarz & Zdziarski 1995) with the reflection fraction 
R — Q/2n left free to vary (where R = 1 denotes reflection 
from material subtending 2n sr as seen from the illuminat- 
ing source). The cut-off energy for the pexrav component is 
fixed at 300 keV for consistency with the built-in cut-off in 
the later used reflionx model (similarly the powerlaw cut- 
off is fixed at 300 keV) and cosi = 0.87 throughout. Narrow 
Gaussians of fixed width (10 eV) are added where appro- 
priate to model neutral and ionized emission or absorption 
features. FeK/3 emission at 7.056 keV is accounted for with 
fixed line energy and flux Fk/s = 0.13 x Fko. The flux of 
the FeK/3 component measured here is frozen and carried 
forward into later models since it is not included within re- 
flionx. Any broad residuals in the FeK region are mod- 
elled using a broad Gaussian with a width free to vary. The 
parameters obtained during this initial parametrisation are 
hsted in Table [H] 

The baseline model is intended to model the entire 0.6- 
100.0 keV range in full to account for features such as any 
soft excess which may be present and distant reflection in a 
more self consistent manner than in the initial parametriza- 
tion. For example, by replacing pexrav with reflionx to 
model reflection and using high ionization xstar grids to 
model absorption in the FeK region. This baseline model 
does not, however, include any components which may be 
representative of emission occurring from regions close to 
the central black hole i.e. any broadened lines or observed 
red-wing are left unmodelled for the purposes of this first 
step, or null hypothesis. 

A Comptonization origin of the soft excess is assumed 
in this model, using COMpTT with a soft photon input tem- 
perature fixed at 0.02 keV. This is quantitively similar to 
adding a steeper soft powerlaw component. The unblurred 
distant reflection component is accounted for initially by 
PEXRAV (plus Gaussians to model the Fe K region) and then 
as described in Section 3.1.2 above for the baseline model. 
Additional soft X-ray emission lines (such as O VII, O VIII 
etc) are added on an ad-hoc basis where required, again with 
fixed width a = 0.01 keV. The warm absorber component is 
treated as stated in Section 3.1.3 with successive zones of an 
XSTAR grid added as statistically required. As expected from 
previous studies of these objects, MCG-6-30-15, NGC 3516, 
NGC 3783 and NGC 4051 all require at least one ionized 
absorption zone. In Fairall 9, however, the Galactic line-of- 
sight absorption is sufficient to describe the total absorption 
at soft X-ray energies. The warm absorber in these objects 
may indeed be more complex than modelled here, although 
the main aim is to derive and measure black hole and accre- 
tion disc properties from the Fe K region, given a consistent 
model to the broadband Suzaku data. Once the baseline 
model has been parametrised we can then proceed to exam- 
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Figure 1. Ratio plots of tlie 4-91teV residuals without any modelling of the reflection component i.e. the FeK region is left totally 
unmodelled with only warm absorption at soft X-ray energies being taken into account in addition to the intrinsic powerlaw and 
COMpTT (where required). The vertical dashed lines represent 6.4 keV in the observed frames. 



ine the nature of the residuals in the FeK region and the 
strength of any broad red-wing indicative of emission from 
the inner regions of the accretion disc. 

3.2.1 Fairall 9 

Fairall 9 is well modelled using an approach similar to that 
used in Patrick et al. (2011), consisting of COMpTT to model 
the soft excess plus powerlaw and an unblurred reflionx to 
account for distant reflection. The broadband spectral prop- 
erties in this long observation are very similar to those in the 
shorter 2007 observation, however, the initial parametrisa- 
tion measures a reflection fraction of 7i = 1.551q'24 com- 
pared to 7? = 0.52l°;^^ in the previous 2007 observation. A 
broad Gaussian improves the flt to the FeK region with 
Ax^ ~ —13 and parameters which are consistent with 
Patrick et al. (2011), see TableED Fitting both the 2010 and 
2007 observations simultaneously and allowing the reflection 



fraction and intrinsic powerlaw to vary gives consistent re- 
sults with both analyses and a smaller difference between the 
reflection fractions for each observation. Forcing the power- 
law F to be the same in both observations provides a good 
fit to the XIS spectra of both observations with the majority 
of the residuals remaining in the fit to the HXD/PIN data. 

Due to the increased signal-to-noise in this observation, 
an additional narrow emission line due to O VIII Lya is ob- 
served at 0.66lo Q2 keV with an improvement of Ax^ ~ —35 
for a further two free parameters. However, since the re- 
flionx model includes some soft X-ray emission lines, in 
this high quality data it is these lines which appear to de- 
termine the obtained refiection parameters and are there- 
fore detrimental to the overall fit. In particular the hard 
X-ray data from the HXD is poorly fit. In order to resolve 
this issue, a neutral absorber at the redshift of the source 
{z — 0.047) is included to only absorb the soft X-ray emis- 
sion lines included within reflionx. This results in an im- 
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provement of Ax^ = —109 in the overall fit with column 
density A^h = 39.59t^;i9 x 10^^ cm"^. Prior to including 
any diskline models indicative of relativistic emission, ion- 
ized narrow emission lines which may result from FeXXV 
and FeXXVl are found at B-GSlf^'o? keV and 7.001^5 keV 
respectively (see Table IX2| . These hues were also found in 
the shorter 2007 observation of Fairall 9 (SchmoU et al. 2009; 
Patrick et al. 2011) with consistent energy and EW. 



3.2.2 MGG-6-30-15 

In this analysis we take the time averaged spectrum of 
MCG-6-30-15 during the January 2006 observations, prehm- 
inary analysis reveals that the dervied FeK region param- 
eters are consistent for all three observations and as such 
they are co-added as per Miniutti et al. (2007). Most ap- 
parent is the presence of a small soft and a hard excess, 
in addition to the requirement for some degree of warm 
absorption. A very poor fit to the data is found without 
any warm absorption with > 60. A single low column 
(A^H ~ 7.5 X 102^cm-^ log^ ~ 1.6) warm absorption 
zone improves the fit to xS ~ 3.2. A second low column 
(A''h ~ 2.4 X 10^^ cm~^) ionized absorption zone improves 
the fit by Ax^ ~ —1284 to x^ ~ 2.5 with both zones having 
ionizations of log^ ~ 0.96 and log^ ~ 1.94. 

In order to parametrise the reflection in MCG-6-30- 
15, we initially use the pexrav neutral refiection model. 
This yields a refiection fraction of i? = 0.79t(! og with a 
good fit to the HXD and BAT data. This then allows us 
to measure the strength of features in the complex FeK 
region prior to employing a more self-consistent refiection 
model such as reflionx. There is strong neutral Fe Ka and 
what may be ionized FeXXV emission improving the fit by 
Ax^ ~ —331 and -51 respectively (see Table IXTj) : it is, how- 
ever, possible the ionized feature may be the blue-wing to 
a broadened FeK line profile. We also find ionized FeXXVI 
emission which is fixed at 6.97 keV throughout the fits. In 
addition to this, very clear in the FeK region are two ion- 
ized absorption features at 6.76to'o2 keV and 7.05to o2 keV 
(see Figure [T]) indicative of the presence of a highly ionized 
absorber (as also found in a Chandra HETG observation 
of MCG-6-30-15 by Young et al. 2005). Subsequent to mod- 
elling these narrow features there is still a poor fit to the data 
(xS ~ 1.3), however there are strong residuals remaining ex- 
tending from ~ 5.0 — 6.4keV. Including a broad Gaussian to 
model this feature results in a line centred at 5.93^0 54 keV 
and EW = 149tg^ eV. Previous analyses of this object at 
energies > 3 keV found a much stronger broadened iron line, 
for example Miniutti et al. (2007) (when modelling with 
a 'double Gaussian' to account for both the red and blue- 
wings) find consistent line energy and (Twidth but a much 
higher cumulative EW = 3201*5 eV in an analysis of the 
same time averaged data used here. 

Replacing the pexrav model and narrow 6.4 keV Gaus- 
sian with REFLIONX (FeK/3 emission with fiux fixed as in 
Table to form the baseline model provides a poor fit to 
the data (x^ ~ 1.82). Despite the two absorption zones and 
soft emission included with reflionx, further soft narrow 
emission and absorption features are required. These can be 
attributed to O VIII Lya, O VII RRC, Ne IX resonance emis- 
sion and SXV/ XVI ls-2p absorption fines at 2.37lo;o} keV 



and 2.77^0 08 keV (see Table IXSl for a summary of narrow 
soft X-ray lines in each object). 

After using reflionx to model the refiection, there is 
still an excess at hard X-ray energies. Since the baseline 
model does not assume any emission from the inner regions 
of the accretion disc, the additional hard excess is therefore 
accounted for using a high column density partial covering 
model to supplement the reflionx hard X-ray component. 
The resulting model provides a reasonable fit to the data at 
both soft and hard X-ray energies with x? ~ 1-26 and an ob- 
scuring column density A^h ~ 3.4 x lO'^* cm~^ resulting in a 
covering fraction of Cfrac ~ 50%. It should be noted that due 
to the high column of the partial coverer, there is some effect 
to the residuals in the FeK region and as such using such 
a model to account for the hard excess infiuences findings 
relating to the strength of any relativistic emission which 
may be present (see Figure [2] and Section 3.3.1 for more de- 
tails). A good fit is obtained to the refiection continuum and 
neutral FeKa line with purely Solar iron abundance. 

The remaining residuals in the Fe K region indicate the 
possibility of a moderate red-wing below 6.4 keV extending 
down to ~ 5keV in addition to narrow ionized emission 
from FeXXVI and strong absorption due to FeXXV and 
FeXXVI. Following the procedure stated in Section 3.1.4, 
we use an xstar generated grid to account for the highly 
ionized absorption, as a consequence of this some of the ob- 
served red-wing may become weaker as noted in Reeves et al. 
(2004). This improves the fit further by Ax^ = -298 for an 
additional 3 free parameters and the residuals show no fur- 
ther absorption. The highly ionized absorption parameters 
found here are consistent with those found in past studies 
of MCG-6-30-15 (Miller et al. 2008; 2009). The inclusion of 
weak FeXXVI Lya emission with a reduction of Ax^ = — 7 
finalises the baseline model. 

The baseline model (Table IA2|) . therefore, provides a 
reasonable fit to the time averaged January 2006 without 
the need for any strong emission from the inner regions of 
the accretion disc with a final fit statistic of xl ~ 1-11. 
However, some previous analyses of this object have found 
no such feature i.e. partial covering and absorption com- 
ponents alone have been claimed to give an acceptable fit. 
For example, Miller et al. (2009) model MCG-6-30-15 with 
'clumpy' absorption zones and achieve a good fit to the data 
without the need for the addition of a diskline component 
representative of emission from the inner regions of the ac- 
cretion disc. In the present analysis, however, we do find that 
there is still some indication of a moderate red-wing which 
may be due to a relativistically broadened component (see 
Figure [2] this is later fit with a kerrdisk line profile with 
an improvement of Ax^ ~ —59, see Section 3.3 and Table 

3.2.3 NGC 3516 

Taking a similar approach as above, we again consider the 
time averaged data from the October 2009 deep observation 
(~ 251 ks) of NGC 3516. Initially modelling the broadband 
spectrum with a simple Galactically absorbed powerlaw and 
a single warm absorption zone, a poor fit is obtained with 
X^ ^ 6.6 and photon index of F ~ 1.7. The most striking 
residuals are at hard energies and the narrow FeKa core. 
There is no indication of a soft excess as found in an analysis 
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Figure 2. Ratio plots of the XIS spectra to the basehne model (i.e. Table lA2l without including a broad FeKa) revealing excesses at 
energies red-ward of 6.4 keV in some objects. Plot MCG-6-30-15 A shows the residuals to the FeK region without the use of a high 
column partial covering component, plot MCG-6-30-15 B shows the residuals of the complete baseline model as per the other objects. 
Note that the high column partial coverer significantly reduces the strength of the red-wing in MCG-6-30-15. 



of the 2005 observation by Markowitz et al. (2008) although 
NGC 3516 is highly variable. 

Parametrising the reflection continuum with PEXRAV re- 
sults in a reflection fraction oi R = 0.87^0 22 ^-nd ^ relatively 
simple FeK region with a very strong narrow FeKa line 
with EW = 229tg^ eV improving the fit by Ax^ = -1396 
in addition to narrow FeXXV emission (Figure [T]). A broad 
Gaussian is added to the model which improves the fit by 
Ax^ 17 and EW = 55^11 eV as noted in Table IXT] 

Replacing the refiection continuum with that from re- 
FLIONX (and no longer modelling the Fe K region with Gaus- 
sians), however, leaves additional excesses at hard X-ray en- 
ergies and requires the introduction of a high column partial 
covering component, consistent with the method used above 
for MCG-6-30-15. This makes a further significant improve- 
ment to the fit (Ax^ ~ —94) with a covering column density 
of A*'h ~ 1.2 X lO'^^ cm~^, ionization log^ ~ 0.56 and a cov- 



ering fraction of Cfrac ~ 21%, resulting in a reasonable fit 
to the broadband data. 

There are residuals at soft energies possibly originating 
from O VII recombination, NelX, MgXII and what may be 
MgXI ls-3p emission (see Table 1X3)) . in addition to these 
there is a narrow emission feature observed at 7.42 j^g gg keV 
due to NiKa with EW — 37IJ4 eV. Analysis of any remain- 
ing residuals in the FeK region do not indicate any highly 
ionized absorption at FeXXV/XXVI energies. Narrow emis- 
sion keV is observed and modelled with a narrow 
Gaussian, improving the fit by Ax^ = —8 for an additional 
two free parameters. This baseline model is a very good fit 
to the data with Xv ~ 1-053 without any relativistic line 
profile. There is little indication of any strong red-wing be- 
low 6.4 keV after full modelling of the broadband continuum 
(Figure [2] & Figure |3]|. However, it is possible that the par- 
tial covering component may be reducing observed spectral 
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curvature which may in fact result from relativistic emission 
(see Section 3.3.1). In a previous analysis a 2005 Suzaku ob- 
servation of NGC 3516, a relativistically broadened line was 
required (Markowitz et al. 2008), however the 2-10 keV flux 
is higher by a factor ~ 1.8 in the 2005 observation. 

3.2.4 NGC 3783 

In NGC 3783 there is a clear soft and hard X-ray excess 
along with some degree of intrinsic absorption, these are 
initially modelled as above using the COMpTT, pexrav and 
warm absorber components in addition to the intrinsic pow- 
erlaw and Galactic absorption (using WABS). Two absorption 
zones are required, including a low column, low ionization 
zone (A^H ~ 1 X 10^1 cm-2, log^ ~ 0.7, see Table R2l) con- 
sistent with Reeves et al. (2004). This still yields a poor fit 
to the data with xi ~ 2.4, however, there are clear residuals 
at soft energies due to narrow emission. 

The observed soft X-ray emission lines are included us- 
ing narrow Gaussians due to O VIIILya, OVII recombina- 
tion, NelX recombination and MgXIILya (see Table 1X3)) . 
In addition to this there is a narrow emission feature at 
3.281q o2 keV which may be due to Ar XVIII, improving the 
fit by Ax^ = -26 and a significant narrow absorption fea- 
ture at 2.39^0 01 keV due to ls-2p SXV with an improve- 
ment of Ax^ = —163 (despite the inclusion of two warm 
absorption zones, however, parameters associated with the 
SXV absorption may be influenced by AuM-edge calibra- 
tion uncertainties). 

Following on from the Reeves et al. (2004) analysis of 
a long XMM-Newton observation of NGC 3783 in 2001, we 
can confirm here the presence of strong FeXXV absorption 
in the FeK region in addition to weaker absorption due to 
less ionized Fe at ~ 6.5 keV (Figure [1]). Modelling this ab- 
sorption with two simple Gaussians results in an improve- 
ment by Ax^ = —38 and Ax^ = —7 respectively. However, 
to form a physical model to account for the observed com- 
plexity in the FeK region we use a high ionization xstar 
grid (as in MCG-6-30-15), whilst not as good a fit as the 
two simple Gaussians, the fit is improved by Ax^ ~ —18. 
The high ionization zone successfully models the FeXXV 
absorption with the parameters given in Table E2l although 
the less ionized Fe absorption at ~ 6.5 keV is left relatively 
unmodelled. Attempting to model this feature with a sec- 
ond highly ionized absorption zone is unsuccessful due to its 
relatively low statistical weight and as such the absorption 
feature at ~ 6.5 keV is left unmodelled. 

An additional narrow Fe XXVI emission line is included 
(Ax^ = -45) to give the baseline model for NGC 3783 
which provides a good fit to the data (xt = 1.083). How- 
ever, the inclusion of a highly ionized absorption zone still 
leaves a moderate red-wing below 6.4 keV extending down 
to ~ 5keV (Figure [2|. This indicates the possible presence 
of statistically significant relativistic emission from the inner 
regions of the accretion disc (later fit with a kerrdisk line 
profile with an improvement of Ax^ = — 32, see Section 3.3 
and Table El). 

3.2.5 NGC 4051 

The X-ray spectrum of NGC 4051 is very complex as well 
as variable, differing between all three observations to date 



with Suzaku. An in depth study of these observations with 
the aid of HETG data from Chandra is conducted in Lobban 
et al. (2011) and we take note of their findings in our mod- 
elling of the 2008 deep observation (Obs ID:703023010) dur- 
ing a period of relatively high flux for this source -F2-10 kov = 
2.45 X 10~^^ergcm~^ s~^. There are clear soft and hard ex- 
cesses which are modelled with COMpTT and pexrav re- 
spectively, along with two low column A*'h ~ 1 x 10^^ cm~^ 
warm absorber zones of the xstar grid with ionization pa- 
rameters log^ ~ 3.0 and log^ ~ 1.8, providing a still poor 
flt to the data {xl ~ 1.5). 

Given the above model of the broad-band data, the re- 
maining residuals lie in the FeK region as can be seen in 
Figure [1] As documented in Pounds et al. (2004) and Lob- 
ban et al. (2011) there is evidence of a highly ionized outflow 
resulting in blue-shifted FeXXV and FeXVVI absorption 
lines at 6.79lQ Q5keV and 7.12lQ g3 keV. A parametrisation 
of these features with narrow Gaussians yields a good fit to 
the data with xt ~ l-li see Table {KT\ 

Replacing the pexrav with the reflionx reflection 
model results in a poor fit to the reflection spectrum, pre- 
dominantly in the HXD and BAT. In order to supplement 
the refiection continuum, we add a high column partial cov- 
erer (Table IA2I) which improves the fit by Ax^ = —62 
(as above in MCG-6-30-15 and NGC 3516, no such ex- 
cesses remain when using pexrav, 7? = 1.52lg'29), result- 
ing in a good fit to both the high energy data and the re- 
flection continuum. We flnd a sub-solar iron abundance of 
Zfc = 0.5lg'4 consistent with Lobban et al. (2011). We also 
flnd an emission feature at 2.22tQ g3 keV with an improve- 
ment of Ax^ = — 16 which may be due to ls-3p Si XIII. 

Lobban et al. (2011) find that the high ionization ab- 
sorption lines in the FeK region are best fit by a highly 
ionized XSTAR grid with a turbulent velocity of 3000 kms^^, 
however, in this analysis we continue to use a turbulent ve- 
locity of only 1000 km s~^ to remain consistent with the pro- 
cedure used above and outlined in Section 3.1.4. This still 
models the absorption features well, improving the fit by a 
further Ax^ = —77 with parameters consistent with Lob- 
ban et al. (2011). Subsequent to modelling the highly ion- 
ized zone and narrow FeXXV emission, there are very few 
residuals remaining below 6.4 keV which indicates that rela- 
tivistic emission from the inner regions of the accretion disc 
may not be required in NGC 4051 (see Figure dj. However, 
we also test for the possibility of a relativistic component 
without the use of a partial covering component in Section 
3.3.1. 

3.3 Kerrdisk profile 

Given the above baseline model in Table [X2I for each object, 
we can accurately examine any residuals in the Fe K region. 
As a step towards assessing the strength of any contribution 
from relativistically blurred emission from the inner regions 
of the accretion disc, we employ the kerrdisk line profile 
model with BH spin as a variable parameter (Brenneman 
& Reynolds 2006). Of course, this is under the assumption 
that the broad residuals remaining in the FeK region are 
due to emission from the inner regions of the accretion disc. 

This allows us to measure the strength of the relativis- 
tic emission, the inner radius of the disc, the inclination 
and the emissivity of the disc. The emissivity index scales 
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Figure 3. uFu plots of the baseline model including any absorption components: powerlaw (green), COMPTT (light blue), REFLIONX 
(purple), partial coverer (yellow), narrow Gaussians to model additional FeK emission (dark purple) and soft emission lines (orange) 
where required. HXD data is in red and BAT data from Swift is represented by blue circles. 



as i?"' where g ~ 3 would be expected from the inner re- 
gions of an accretion disc and q > 5 indicates that emission 
from the accretion disc is very centrally concentrated with 
a significant amount of hght bending. We assume the ac- 
cretion disc extends down to the inner-most stable circular 
orbit (-Risco), which itself varies according to the value of 
the BH spin parameter. We fix i?out = 400 Rms and assume 
relativistically blurred neutral Fe Ka with rest frame energy 
6.4 keV. The kerrdisk model only considers the possibility 
of prograde BH spin with < a < 0.998 where a = cJ/GM^ 
(J is angular momentum). 

It should be noted that the results obtained here can 
only be judged as estimates due to the possible model degen- 
eracies. In the event that residuals in the Fe K region result 
from the observation of relativistically blurred diskline emis- 
sion, we attempt to consider these degeneracies and present 
our best estimates of the line parameters. When employing 
a diskline model (e.g. kerrdisk), care should be taken due 



to the interplay between the emissivity index q and the spin 
parameter a. As such the spin versus Ax^ plots in Figure 
[4] indicate the confidence levels for two interesting parame- 
ters where both q and a are allowed to vary. The presence 
of narrow ionized emission lines in the Fe K region can also 
affect the measured emissivity index q and hence the value 
of the BH spin obtained with a diskline model. However, in 
these deep observations the number of counts is generally 
sufficient to distinguish between the presence of a narrow 
Gaussian or the blue-wing of a relativistically blurred line 
profile. 

The analysis of the long 2010 Suzaiu Fairall 9 observa- 
tion here yields, similar results to those of the shorter 2007 
observation in Patrick et al. (2011). The disc emissivity is 
consistent with q = S.Slg 4 although despite the higher sig- 
nal to noise data, it is not as well constrained (compared 
to q = 2.7lg'4). analysis by SchmoU et al. of the 

same 2007 data they place only a lower limit of g > 4.9. 
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Table 3. kerrdisk model parameters in addition to the baseline model. Absorption and OOMPTT components are consistent with 
those quoted in Table IA2I and are not stated again here to avoid repetition. * Denotes a frozen parameter. " Flux given in units 
(10~^ phcm~^ s~^). * Denotes a frozen parameter. A positive Ax^ represents a worsening in the fit upon the removal of the component. 
Spin cannot be constrained in NGC 3516 at the 90% confidence level. 
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A significant difference which is responsible for the poorly 
constrained emissivity index is the clearly detected narrow 
Fe XXV emission which is now included in the fit to the Fe K 
region. The interplay between the narrow Gaussian used to 
model the feature at ~ 6.7 keV and the blue- wing of the ker- 
rdisk line profile causes degeneracy between the emissivity 
index and FeXXV flux. The inclusion of narrow FeXXV 



emission also results in a change to the inferred value of the 
spin parameter a, here it is estimated to be a = 0.67lo }i> 
whereas a lower value of a = 0.441q5i is obtained in the pre- 
vious analysis when the narrow emission was not included 
(Figure |4]). The flndings made here are, however, consistent 
with an analysis of a 130 ks XMM-Newton observation of 
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Figure 5. Ratio plots of the basclinc+KERRDlSK model between 4-9 keV. After the application of the kerrdisk relativistic line emission 
model there are few remaining residuals below 6.4 keV, giving a very good fit to the Fe K region and the broad-band data as a whole. 



Fairall 9 in 2009 by Emmanoulopoulos et al. (2011) who 
find a = Q.iQtl fo- 

Maximal spin is ruled out in all objects, including 
MCG-6-30-15 in contrast to claims to the contrary by 
Reynolds et al. (2005) and Miniutti et al. (2007) of a > 0.93 
and a > 0.92 respectively. One possible difference is that 
a full broad-band model has been compiled to fit the 0.6- 
100.0 keV spectrum, whereas the Miniutti et al. (2007) spin 
measurement is only valid over the 3-45 keV range and with- 
out including warm absorbers. In this analysis we find a 
fairly typical emissivity index of g = 2.7lo^, firstly sug- 
gesting that emission in MCG-6-30-15 is not particularly 
centrally concentrated and the effects of light bending may 
be minimal. The spin of the central black hole in MCG-6- 
30-15 is estimated to a much lower spin a = 

0-49lo:?2 (Table 

[3ll, maximal spin is ruled out at > 99.5% confidence in this 
model where the fit is made worse by Ax^ = +70 for the 
case when the spin parameter is fixed at a = 0.998 (see 
Figure |4]). 



Spin cannot be constrained in NGC 3516 and only loose 
constraints can be made on the emissivity index of the disc 
q = 2.5I1Q, this is expected due to the relatively weak line 
modelled using a broad Gaussian in Table lAll The fit is 
only improved by Ax^ = — 10 after adding kerrdisk emis- 
sion with EW — 46l3g and no further narrow ionized 
emission is required in the FeK region with the kerrdisk 
profile taking the place of the narrow Fe XXV emission (see 
Figure [SJ. However, the FeK region is well described by the 
baseline model (Table IA2p which includes narrow Fe XXV 
emission, adding a kerrdisk line profile only improves the 
fit in this case by = —2 over the baseline model. 

Only an upper limit of a < 0.32 at the 90% confidence 
level can be placed upon the spin in NGC 3783 using the 
KERRDISK model. A maximally rotating BH is again ruled 
out with a worsening to the fit of Ax'^ ~ +15. Similarily 
only an upper limit can be placed upon the inclination of 
the accretion disc with i < 17°, however, the emissivity can 
be tightly constrained to q = 2.61q'^. The strength of the 
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broad line component is consistent with Reeves et al. (2004) 
{EW = 46^3^ eV here compared to EW = SSl^a eV), being 
reduced significantly after including a physical model for the 
highly ionized absorption. 



The broad line feature in NGC 4051 makes a statisti- 
cally insignificant improvement to the quality of the fit with 
a reduction of Ax^ = —3. Hence, none of the accretion disc 
parameters can be particularly well constrained, the emis- 
sivity index is fixed at g = 3 and only an upper limit can 
be made upon the spin parameter of a < 0.94. The FeXXV 
emission listed in Table IX2l is no longer required with the in- 
troduction of the KERRDISK component and is alternatively 
explained as a blue- wing to the line profile (Figure |6} . The 
results obtained here using a highly ionized xstar grid with 
a turbulent velocity of 1000 kms~^ are consistent with those 
when a turbulent velocity of 3000 km s"'^ is used (e.g. by 
Lobban et al. 2011) i.e. line strength, EW and _Rin. 



3.3.1 The effect of partial covenng 

As stated in Section 3.2.2, the high column density partial 
covering component affects predominantly the hard X-ray 
data in the HXD and BAT detectors, however, we test the 
effect of this partial covering component on any curvature in 
the Fe K region. The partial covering component and distant 
reflection from reflionx are removed and replaced with 
the PEXRAV neutral reflection model (plus a narrow neu- 
tral 6.4 keV Gaussian) which achieved a good fit to the data 
in the initial parametrisation. This interpretation i.e. with- 
out partial covering, yields a slightly higher spin constraint 
of a = 0.82l^;J^g for MCG-6-30-15, however, a maximally 
rotating BH is still discounted at an increase of Ax^ ^+17 
(> 99.5% confidence) despite a fairly low refiection fraction 
Figure [7] shows the residuals to the data 
with the partial covering component removed (still including 
reflionx), showing that the high column partial covering 
component in MCG-6-30-15 predominantly affects the fit to 
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Table 4. kerrdisk model parameters when the reflection con- 
tinuum is modelled with PEXRAV in place of the REFLIONX plus 
partial coverer configuration i.e. no partial covering component is 
used here. " Flux given in units (10~^ phcm^^ s^^). * Denotes a 
frozen parameter. A positive Ax^ represents a worsening in the 
fit. Spin cannot be constrained in NGC 3516 at the 90% confi- 
dence level. 
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the hard X-ray data. Only subtle curvature is introduced to 
the 7-10 keV spectrum in the XIS detectors. 

Both NGC 3516 and NGC 4051 (which display no 
strong broadened emission) also feature a partial covering 
component in the baseline model, removing the partial cov- 
erer and instead modelling the hard excess using pexrav 
gives accretion disc parameters which are entirely consis- 
tent with those in Table [3] (see Table |4]). This again suggests 
that significant relativistically broadened emission is prob- 
ably not strongly present in these objects and the high col- 
umn partial coverer does not have a significant affect upon 




10 20 
Observed Energy (keV) 

Figure 7. Comparison between the fit to the 4-60 keV data for 
MCG-6-30-15 with and without a partial covering component. 
The first plot shows the baseline model, in the second plot the 
partial coverer has been removed and the data refit i.e. still using 
REFLIONX and not PEXRAV. Note that the (high column) partial 
covering component therefore affects predominantly the hard X- 
ray data, with little affect upon the Fe K region. HXD data is in 
red and BAT data from Swift is represented by blue circles. 



the Fe K modelling. The partial covering component as used 
in the baseline model does have some affect upon the derived 
accretion disc and BH parameters, although only in MCG-6- 
30-15, for example the measured strength (i.e. flux and EW) 
is slightly higher without the presence of a partial coverer 
as is the measured BH spin (see Table HJ. Properties of the 
accretion disc such as the emissivity index and the inclina- 
tion remain relatively unchanged. Despite the removal of the 
partial covering component in all objects (where present), a 
maximally spinning central BH is ruled out in these objects. 
Note that this is not to say that the residuals in the FeK 
region cannot be accounted for by a partial coverer, sim- 
ply that the high column partial coverer used here does not 
significantly affect the measured accretion disc parameters. 



3.,S.2 Testing for retrograde SMBH spin 

Recent work has suggested the possibility of retrograde 
SMBH spin in AGN, whereby the rotation of the central BH 
and the surrounding accretion disc are anti-parallel whilst 
maintaining a stable configuration (King et al. 2005). The 
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RELLINE model (Dauser et al. 2010) is a relativistic line 
emission model accomodating both negatively and positively 
spinning central BHs, which produces line profiles consis- 
tent with established prograde spin line emission models 
e.g. KERRDISK (Brenneman & Reynolds 2006) and kyrline 
(Dovciak et al. 2004). 

The spin versus Ax^ plot for NGC 3783 in Figure [4] 
reveals a best-fitting spin parameter of a = according the 
KERRDISK model, assuming a strictly prograde BH spin, in- 
deed NGC 3783 is the only object showing this trend in this 
sample. We use the relline model in addition to the base- 
line model to assess the suitability of a retrograde BH in a 
fit to the FeK region in NGC 3783. In a similar procedure 
to that employed using the kerrdisk model, we assume the 
accretion disc extends down to -Risco and that emission oc- 
curs from neutral Fe Ka and as such we fix the line energy 
at 6.4 keV. 

The accretion disc parameters obtained are very similar 
to those obtained with kerrdisk in Table |3] with a slightly 
higher emissivity oi q = 3.21q4, inclination i < 17° and 
EW = 43lfgeV (see Table [S]). Modelling the FeK region 
of NGC 3783 allowing for retrograde spin suggests an upper 
limit of a < —0.04, as shown in Figure [S] This modelling of 
NGC 3783, therefore, suggests a negatively spinning BH in 
NGC 3783 at the 90% confidence level, ruling out a max- 
imally spinning prograde BH at > 99.5% confidence with 
an increase of Ax^ ~ -fl3 in fit statistic when fixed at 
a = -1-0.998. Figure [8] also shows a 2D emissivity vs spin 
contour plot for NGC 3783 using relline, indicating the 
relation between these two parameters. It is interesting to 
note that the emissivity of the accretion disc is well con- 
strained and, in this object, has little effect upon the mea- 
sured BH spin. There is only a slight negative correlation 
in which emission is more centrally concentrated for lower 
(more negative) BH spin. 

It should be noted that while a retrograde BH does 
indeed push the innermost stable circular orbit of the ac- 
cretion disc outwards, thereby reducing the extent to which 
the line profile can become relativistically broadened, the 
line profile is noticibly broadened and not simply represen- 
tative of a series of narrow Gaussians which would otherwise 
be indicative of purely distant emission (see Figure O . 

3.4 Dual reflector 

The previous kerrdisk relativistic line emission model is 
representative of refiection off the inner regions of the accre- 
tion disc. If strong emission arises from this region, the same 
relativistic broadening must also be applied to the reflec- 
tion continuum as a whole, as opposed to purely the narrow 
6.4 keV FeKa emission. To account for this we replace the 
kerrdisk model with a second reflionx component con- 
volved with the KERRCONV kernel (Brenneman & Reynolds 
2006) in addition to the unblurred reflection component. We 
also remove the partial covering component from the dual 
reflector fits to assess this alternative interpretation. 

In some analyses of AGN, such a blurred reflection com- 
ponent has been used to model the soft excess under the 
assumption that the regions responsible for the smearing of 
discrete atomic emission features at soft energies are also 
responsible for broadening in the FeK region (e.g. SchmoU 
et al. 2009; Nardini et al. 2010). However, in an albeit 




° -0.5 0.5 

Spin Parameter a 

Figure 8. The upper plot is a fit statistic versus a plot for SMBH 
spin from the relline model when used in conjunction with the 
baseline model in NGC 3783. The spin parameter a is allowed to 
cover the full range —0.988 < a < +0.998 to assess the possi- 
bility of retrograde BH spin. The lower plot shows a 2D contour 
plot showing the relation between emissivity index and spin, con- 
fidence levels labelled in both plots are 68% (black), 90% (red), 
99% (green) and 99.9% (blue) for two interesting parameters. 



small sample of bare Seyfert AGN, Patrick et al. (2011) de- 
rived contrasting accretion disc parameters depending upon 
whether such blurring is allowed to model the soft excess 
in addition to features in the FeK region. A signiflcant 
amount of blurring is required in order to smear the soft 
atomic emission features into a smooth soft excess, resulting 
in much higher derived BH spin and very centrally concen- 
trated emission (i.e. high g > 4). 

Taking this into account, in this interpretation of the 
data, we retain the use of the COMpTT Comptonization ori- 
gin to model the soft excess, leaving the convolved refiection 
spectrum to account for any observed iron line broadening 
and hard excess over purely distant reflection. A purely re- 
flection origin for the hard excess is assumed and no partial 
covering model is used in the dual reflector flt as in MCG- 
6-30-15, NGC 3516 and NGC 4051. The iron abundance of 
both the inner and distant reflector are tied and frozen at 
the best fltting values from the baseline model through the 
flts. 

A good flt can be obtained to all objects with a dual re- 
flector model (Table EH), with the exception of MCG-6-30- 
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15, the hard X-ray data from the HXD and BAT detectors 
are fit well without the requirement for any additional reflec- 
tion or a partial covering component. Due to the increased 
number of free model parameters and hence complexity in 
the broadband model it is difficult to measure accretion disc 
parameters (particularly BH spin) to any reasonable accu- 
racy, even with high quality data as used here from deep 
Suzaku observations. Consequently, no better constraints 
can be made using a dual reflector fit, rather the results ob- 
tained are generally consistent with the baseline-f kerrdisk 
model (Table [3|. 

For the case of MCG-6-30-15, a reasonable quality of fit 
{xl = 2061.4/1823) is obtained with a dual reflector fit. As 
noted above, the interpretation of the soft excess can dra- 
matically alter the measured accretion disc parameters. In 
an object as complex as MCG-6-30-15 with strong reflection, 
broadened relativistic emission in the FeK region, complex 
warm absorption and a soft excess, the use of a dual reflec- 
tor results in a number of degeneracies when modelling the 
spectrum. As a result of this we fix the COMpTT parameters 
at the best-fit values obtained with the baseline model and 
then proceed to add a blurred reflection component in place 
of the partial covering model. This approach ensures that 
the blurred reflection component does not model the soft 
excess in place of COMpTT and as such the accretion disc 
emissivity and BH spin will not be forced to extreme values 
in order to reproduce a smooth and featureless continuum. 

The measured parameters for MCG-6-30-15 are consis- 
tent with those obtained in Table [3] with q = 2.31q j a^nd 
a = 0.611q Also requiring a FeXXV emission line which 
improves the flt by Ax^ ~ —7, however it is weak with a low 
EW — 6^3 eV. Solar abundances are assumed throughout 
and a consistent ionization parameter (^ < 11 erg cm s~^) is 
measured in both the blurred and unblurred reflectors. The 
warm absorber parameters measured with the dual reflector 
fit are also consistent with those obtained in the baseline 
model. 

The properties of the accretion disc are not greatly af- 
fected by the interpretation of the origin of the hard excess 
i.e. whether a high column (A'h ~ 10^'* cm~^) partial cover- 
ing model or a dual blurred and unblurred reflector is used 
(see Table 1X4)) . The partial covering model used in the base- 
line model therefore has a limited effect upon the accretion 
disc parameters and BH spin measured from the Fe K region. 
The results obtained with the dual reflector fit in NGC 3783 
are completely consistent with those measured with the pre- 
vious baseline-|-KERRDlSK model, again an upper limit on the 
BH spin parameter is found at a < 0.31 and the emissivity 
of the accretion disc is also a similar q = 2.81q'2. 

The complex absorption in these objects (other than 
Fairall 9) can mask the extent to which the smeared dis- 
crete emission features are present in the data at soft X-ray 
energies. More extreme blurring (a high emissivity, for ex- 
ample) is not required since the X-ray spectrum at soft en- 
ergies in these more complex objects is far from smooth and 
the discrete emission features do not require smearing into 
a smooth continuum-like shape. As a consequence of this, 
the main components driving the dual reflector fits are the 
features in the Fe K region, rather than curvature due to the 
blurring of soft X-ray atomic emission. 



4 DISCUSSION 

This sample of Seyfert 1 AGN includes deep Suzaku ob- 
servations with exposures > 200 ks at high S/N, enabling 
us to fully model the broad-band spectrum when combined 
with additional hard X-ray data from Swift/BAT. In order 
to robustly assess the nature of any features in the FeK 
region, whether it be absorption or relativistically blurred 
emission, a complete modelling of the X-ray spectrum was 
undertaken to ensure that all features are taken into account. 
With complex X-ray spectra such as found in these AGN, 
deep exposures are essential to help constrain a number of 
features prior to investigating accretion disc properties and 
the potential for measuring BH spin. 

Some improvement is made to the fit for all objects 
with the addition of the kerrdisk model which represents 
blurred line emission from the inner regions of the accretion 
disc. The overall improvement to the fit for NGC 3516 and 
NGC 4051 is smaU with Ax^ = -10 and -3 in each for 4 and 
3 additional parameters respectively, thus even the presence 
of disc emission lines in these AGN is very tentative. Sta- 
tistically much stronger broad line components are found in 
MCG-6-30-15 and NGC 3783 even after a complete mod- 
elling of the warm absorber components, although in some 
circumstances complex absorption models may be able to 
account for these features. 

Narrow ionized emission is found in three of the objects 
here after including a kerrdisk line profile, however the 
blue-wing of the relativistically broadened line profile can 
easily become degenerate with such emission. For example, 
the emissivity index and the inclination of the accretion disc 
can be forced to values such that a narrow feature which 
would otherwise be attributed to FeXXV/ FeXXVI emis- 
sion could be explained by the blue-wing of relativistically 
blurred diskline emission from neutral FeKa. Indeed, NGC 
3516 and NGC 4051 require no ionized narrow emission fea- 
tures in the final baseline-|-KERRDlSK since the broad feature 
can successfully reproduce both the red-wing and narrow- 
like features in the FeK region, yet there is little require- 
ment for broadened line profiles in either object. Emission 
from FeXXV is required only in Fairall 9 even after the 
possible degeneracy mentioned above, FeXXVI Lya emis- 
sion is more common and is found in 3/5 objects (Fairall 
9, MCG-6-30-15 and NGC 3783). Future calorimeter obser- 
vations with Astro-H (Takahashi et al. 2010) will be able 
to distinguish and resolve narrow ionized emission from the 
blue- wing of a diskline profile. 



4.1 Average parameters of the disc 

Given the measured values from the kerrdisk line profile, 
we can calculate typical parameters of the accretion disc. 
The average inclination of the disc is i = 37° ± 8° , although 
this is a small sample, it is consistent with Nandra et al. 
(2007) who find i = 38° ± 6° . From the baseline-|-KERRDlSK 
model (Table [3]) we also measure an average accretion disc 
emissivity of 5 = 2.8 ± 0.2 which (similarly to Patrick et al. 
2011) is a lower value than is often assumed from light bend- 
ing models (e.g. Gallo et al. 2011; Nardini et al. 2011). The 
initial parametrisation in Table lAll with a broad Gaussian 
measures an average strength of the broad line component of 
EW = 77±23 eV at a centroid line energy of 6.11±0.10 keV, 
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also consistent with the previous findings by Nandra et al. 
(2007) with XMM-Newton. Employing the more physical 
KERRDISK model increases the strength to EW = 99±28eV. 

The moderate emissivity indices found here are likely 
due to the accretion disc parameters being obtained purely 
from broadening in the FeK region rather than requiring 
the inner regions of the accretion disc to blur the reflection 
continuum substantially in order to model the soft excess. 
Here we assume the soft excess originates from the Comp- 
tonization of UV seed photons from the disc in a hot corona 
above the disc via the COMpTT model (Titarchuk 1994), 
a second soft powerlaw or disc blackbody would provide a 
quantitively similar flt to the data. If discrete atomic fea- 
tures present in the reflection continuum at soft X-ray ener- 
gies were allowed to model the model the soft excess through 
relativistic smearing due to originating from regions close to 
the central BH, the emissivity of the accretion disc would be 
forced to much higher values i.e. q > 4.5 (see Patrick et al. 
2011). Such high values are not required for a good flt to the 
FeK region alone, therefore, if a blurred reflection compo- 
nent is responsible for the soft excess it must originate from 
a different part of the disc responsible for broadening in the 
FeK region. 

The dual reflector flt in Section 3.4 in which two sep- 
arate REFLIONX components are used for distant and inner 
reflection is also consistent with the typical accretion disc 
measurements made above (Table 1X4)) . While higher reso- 
lution spectroscopy would be required to fully distinguish 
between the varying scenarios, the hard excess can be ex- 
plained by either a high column partial coverer or a dual re- 
flector, although crucially still without the requirement for 
extreme parameters such as high q and a for very centrally 
concentrated emission with an average q = 2.4 ± 0.2. 

4.2 The effect of absorbers on the FeK region 

It is known that narrow ionized emission lines from Fe XXV 
and FeXXVI can alter properties such as the measured in- 
clination and emissivity of the accretion disc, however, the 
warm absorber is often overlooked. For instance warm ab- 
sorber components can introduce spectral curvature, due to 
the bound- free L/K-shell transitions of abundant elements 
(Miller & Turner 2009). In particular, the warm absorber 
can play a signiflcant role in the shape of the spectrum, not 
only at low X-ray energies, but also higher energies closer to 
the FeK region (Miller et al. 2008, Zycki et al. 2010). Util- 
ising the soft X-ray energies below 3 keV allows a complete 
modelling of the warm absorber, adding the corresponding 
spectral curvature, which may not otherwise be taken into 
account without using the available data in the 0.6-3.0 keV 
range. Given a broad-band model, the true nature of the 
FeK region and any broad red-wing can be determined af- 
ter the reduced spectral curvature and observed 'broaden- 
ing' from the required warm absorbers has been taken into 
account. Inappropriate modelling of the warm absorber can 
thus alter the observed strength and dominance of the broad 
red-wing over the continuum, which in turn will affect the 
extent to which black hole spin is required and its measured 
value. 

In addition to the warm absorber properties, modelling 
the highly ionized absorption lines in the FeK region with 
a physically consistent xstar grid (such as with a turbu- 



lent velocity of 1000 km s"'^ used here) ensures the corre- 
sponding curvature is removed from the spectrum. Highly 
ionized zones of gas are required in three of the objects 
here: MCG-6-30-15, NGC 3783 and NGC 4051, also con- 
sistent with past studies (Reeves et al. 2004; Miller et al. 
2008; Tombesi et al. 2010; Lobban et al. 2011). Absorp- 
tion lines in the FeK region due to highly ionized species 
of iron (FeXXV and/ or FeXXVI) are found at the ap- 
propriate energy (or blue-shifted such as in NGC 4051 and 
to a lesser extent in MCG-6-30-15, with outflow velocities 
of 56OOI700 kms'^ and 32OOI500 kms"^ respectively). NGC 
3783 shows both broadened line emission and absorption 
lines in the FeK region, however when a 1000 kms~^ high 
ionization xstar grid is included, the observed strength and 
shape of the broad line is slightly reduced with kerrdisk 
EW = 67tll eV reducing to EW = 45^11 eV (ahhough stiU 
consistent within errors). The same can be said for the highly 
ionized absorption present in MCG-6-30-15 with diskline 
models proving marginally weaker after the introduction of 
the high ionization XSTAR grid (i.e. EW = 180^54 eV re- 
duced to 16lt4^eV). 

As discussed in Section 3.3.1, complex partial covering 
models can provide an alternative explanation of the broad- 
band spectrum and therefore the residuals in the FeK re- 
gion. For example. Miller et al. (2009) construct a model 
which can explain all the features in the broadband spec- 
trum of MCG-6-30-15 with no requirement for a diskline 
component, rather the apparent 'broad' residuals in the Fe K 
region are accounted for with 'complex clumpy' absorbers. 
Here, however, we flnd that even when we apply partial cov- 
ering absorption models to the X-ray spectra, the addition 
of a broad iron line component still appears to be required 
in the case of MCG-6-30-15, but only marginally in NGC 
3516 and NGC 4051. 

4.3 Is MCG-6-30-15 not maximal? 

Previous analyses of MCG-6-30-15 have claimed a very 
strong, broad iron line whilst requiring that the central black 
hole is spinning very rapidly in order to account for the ex- 
treme width of the broad line. Miniutti et al. (2007) measure 
a > 0.92 in a flt which does not consider data below 3keV 
in order to avoid the need to flt for the warm absorber. Here 
we make use of the full broad-band data available, fltting 
in the 0.6-100.0 keV range, using a combination of the XIS, 
HXD and BAT detectors in order to properly account for 
any absorption components which may be present in addi- 
tion to the reflection hump above 10 keV. Taking the warm 
absorber components into account prior to considering emis- 
sion from central regions close to the BH reveals that while 
the majority of the effects due to the warm absorption oc- 
cur at low energies (< 3keV), signiflcant spectral curvature 
is added in MCG-6-30-15 at higher energies (Zycki et al. 
2010). Much of the very broad red-wing observed by Mini- 
utti et al. (2007) is reduced, resulting in a much weaker 
broad line component than has previously been suggested. 
This results in a broad line feature which may originate from 
emission close to a BH with spin a = 0.491q'^2 and lower 
emissivity of g = 2.71q'i while the near maximal spin re- 
quired by Reynolds et al. (2005) and Miniutti et al. (2007) 
is excluded at high confldence levels in our model (see Table 
[3] and Figure |4ll. 



© 2010 RAS, MNRAS 000, ??-?? 



18 A. Patrick et al. 




Observed Energy (keV) 



Figure 9. Line profiles using relline with rest energy 6.4 keV 
and parameters q = 3.0, i = 40° and spin increasing left to right 
a = +0.998 (black), a = +0.5 (red), a = 0.0 (green), a = -0.5 
(dark blue) and a = —0.998 (light blue). Note that even in the 
case of a maximally rotating retrograde BH (furthest out -Risco)i 
there is still a significantly broadened line profile. 

In this paper, hard X-ray emission in excess over those 
which can be reproduced by a single reflionx are accounted 
for using a high column partial coverer. Indeed in MCG-6- 
30-15 Miller, Turner & Reeves (2008) successfully model the 
entire 0.5-45.0 keV spectrum through the use of absorbers 
and partial covering. The full model we form here also in- 
cludes such a component, although its effect upon the broad- 
ening in the Fe K region is subtle. However, even if no partial 
covering model is included, a maximally rotating BH is still 
ruled out at a high confidence level as can be seen in Table 
[4] and Figure [T] 

4.4 A retrograde SMBH in NGC 3783? 

Retrograde SMBH spin has been thought of as a mechanism 
through which the central BH can power radio jets in AGN. 
Garofalo (2009) notes that accretion onto retrograde black 
holes can provide a powerful mechanism resulting in higher 
Blandford-Znajek luminosities than those in the prograde 
spin case due to 7?isco being located further away from 
the SMBH, indeed some radio-loud AGN have been found 
to have inner disc truncations at ~ 8 — 9 7?g e.g. (3C 120, 
Kataoka et al. 2007) with iiisco ~ 9i?g for a maximally 
spinning retrograde BH. 

The line profile produced from emission surrounding 
a negatively spinning BH is narrower than that obtained 
from a positively spinning BH, however measuring a ret- 
rograde spin is not necessarily representative of emission 
occuring from simply distant regions. The inner radius of 
emission of the accretion disc varies according to the spin 
of the central BH with a maximally rotating prograde BH 
producing i?isco ~ 1.23 i?g, a Schwarzschild BH producing 
J^isco = 6-Rg and maximally rotating retrograde BH pro- 
ducing -Risco ~ 9 iig i.e. /?isco increases with decreasing 
positive spin. Even assuming a maximally rotating retro- 
grade BH (the largest value of J?isco), emission still occurs 
close enough to the central black hole for substantial rela- 
tivistic effects and Doppler broadening. The resultant line 



profile is not therefore equivalent to purely distant emis- 
sion. Figure [9] shows a comparison between the line profiles 
produced with the relline model with spin varying from 
a — —0.998 to a = +0.998. However, it could also be that 
the BH is not retrograde, rather that the inner radius of the 
accretion disc is truncated short of the ISCO. 

In this analysis we find a BH spin in NGC 3783 of 
a < —0.04 at the 90% confidence level (for one interest- 
ing parameter), with Figure [S] suggesting a best-fitting in- 
termediate retrograde BH spin. Accretion disc parameters 
obtained with the relline model are consistent with those 
obtained with the kerrdisk model, measuring an emissivity 
q — 3.2tQ'4 and a low inclination of i < 17°. 

A recent analysis of NGC 3783 by Brenneman et 
al. (2011) of the same Suzaku data has suggested a 
BH spin of a > 0.98 with an accretion disc emissiv- 
ity of g = 5.2^0 8 using a dual reflector model i.e. pex- 
mon+(reflionx*relconv). This contrasts significantly 
with the estimates made here, despite Brenneman et al. 
(2011) including warm absorber components in their fit to 
the 0.7-45.0 keV data (while ignoring 1.5-2.5 keV). Brenne- 
man et al. (2011) also allow for partial covering which is not 
required in the baseline model formed here (Table E2)) . The 
highly ionized absorber in the Fe K region detected here (e.g. 
FeXXV absorption, Ax^ = -33) and by Reeves et al. (2004) 
in XMM-Newton data is not included in their fit, instead a 
warm absorber with log^ ^ 2.53 is claimed to model this fea- 
ture sufficiently. There are also variations in the modelling of 
the accretion disc itself, for example Brenneman et al. (2011) 
require a super-solar iron abundance of Zpo = 3.7lo'9 and 
a broken emissivity law with gin = 5.2lg g, rbr = 5.41q g Rg 
and gout ~ 2.91q'2, whereas the analysis here provides a 
good fit to the broad-band data with purely Solar iron abun- 
dance and a single (unbroken) emissivity index of g = 3.2 j^JJ ^ 
(TableO using relline, xl = 1.057). 

In order to test our conclusions here, we form a 
model using essentially the same reflection components 
i.e. PEXRAV+Ka+K/3+(REFLiONX*RELCONv) modified by 
warm absorption, but maintaining Solar abundances and 
the single emissivity law. This brief test confirms the results 
found here and changes the estimated spin to a = — 0.2lg gg 
(xt = 1.072). Further to this, allowing for a broken emis- 
sivity index and allowing the iron abundance to vary does 
not reproduce the Brenneman et al. (2011) results and re- 
laxes the estimates upon the spin parameter a, while the 
abundance remains approximately Solar and there is no in- 
dication for the requirement of a broken emissivity i.e. no 
significant improvement in the quality of the fit. It is possible 
that the treatment of the absorbing components, soft excess 
and assumptions regarding the super-Solar iron abundance 
and broken emissivity law may lead to the higher spin es- 
timated by Brenneman et al. (2011). The overall fit to the 
0.7-45.0 keV data in Brenneman et al. (2011) is poor, mostly 
influenced by the soft X-ray energies (xt = 1.381, in con- 
trast to the best-flt model here of xt = 1.057). However, 
it is not clear whether the narrow soft X-ray emission lines 
present in the Suzaku data are included within the Brenne- 
man et al. (2011) analysis. Indeed a complete modelling of 
the broad-band data (including any soft emission lines) may 
alter the BH spin estimate. 
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4.5 Summary of Suzaku observations so far 

Based upon the results of this paper and those found from 
Suzaku observations of Seyfert Is with low intrinsic absorp- 
tion in Patrick et al. (2011), we can begin to draw conclu- 
sions regarding the FeK region and the accretion disc and 
BH properties. We compare only these two samples to en- 
sure full and consistent broad-band models are used and 
as such estimates of average parameters can be made. The 
objects we consider are therefore: Ark 120, Fairall 9, MCG- 
2-14-009, MCG-6-30-15, Mrk 335, NGC 3516, NGC 3783, 
NGC 4051, NGC 7469 and SWIFT J2127.4-f5654. Results 
are based upon modelling any relativistic broadening in the 
FeK region using the kerrdisk model i.e. Table [3] here and 
Model D (Table 7) in Patrick et al. (2011). 

We find that a maximally rotating Kerr BH is ruled out 
in all objects other that NGC 3516 in which no spin con- 
straints can be made at all upon the spin for the 2009 Suzaku 
observation. Intermediate spin constraints can be made for 
Fairall 9 (a = 0.67t°i;), MGC-6-30-15 (a = 0A9t°ll), 
Mrk 335 (a = 0.70t° J?), NGC 7469 (a = 0.691;^;°^) and 
SWIFT J2127.4-h5654 (a = 0.70^o:l4)- The remaining five 
AGN in this sample have central BHs with a measured spin 
which is consistent with a non-rotating Schwarzschild BH: 
Ark 120 (a < 0.94), MCG-2-14-009 (a < 0.88), NGC 3516 
(no constraints upon a), NGC 3783 (a < 0.32) and NGC 
4051 (a < 0.94). It is possible that a retrogradely spinning 
BH may be found in some of these objects with an upper 
limit upon a, indeed testing the relline model in an anal- 
ysis of NGC 3783 has suggested this with a < -0.04. This 
possibility will be investigated in a subsequent paper. 

Relativistic emission is statistically formally required 
in all but NGC 3516 and NGC 4051 i.e. in 80% of the ob- 
jects, while such emission in MCG-2-14-009 is only weak 
with Ax^ = —12. Although due to the small size of this 
sample (10 objects) limits the extent to which we can com- 
pare this fraction to larger samples and the Seyfert 1 AGN 
population as a whole, it is encouraging that Nandra et al. 
(2007) find broad lines in ~ 70% of objects. However, de la 
Calle Perez et al. (2010) find only ~ 36%, albeit in a flux 
limited sample. 

Over the sample of the 10 objects listed above, we mea- 
sure average parameters of q — 2.5 ± 0.2, a — 0.65 ± 0.05, 
j = 37° ± 3° and EW = 115 ± 16 eV. These compare 
well to Nandra et al. (2007) who find i = 38° ± 6° and 
EW = 91.3 ± 12.8 eV, while de la Calle Perez et al. (2010) 
measure an average emissivity index of g = 2.4 ± 0.4. Al- 
though the average parameters obtained here are by no 
means representative of a larger sample of AGN, each has 
been modelled in a consistent manner. In future work we 
aim to expand this analysis to all the Seyfert 1 AGN in the 
Suzaku archive with a similar uniform broad-band analysis. 

Other objects in the literature which authors have been 
able place constraints upon BH spin include Mrk 79 (Gallo 
et al. 2011) and RBS 1124 (Miniutti et al. 2010). Rela- 
tivistically broadened emission in the Fe K region, although 
occuring outside 6 7?g, has also been found in: MCG-5- 
23-16 (Reeves et al. 2007), MCG+8-11-11 (Bianchi et al. 
2010), Mrk 509 (Ponti et al. 2009) and NGC 2992 (Yaqoob 
et al. 2010); with no broad line component being found 
in: 1H0419-577 (Turner et al. 2009), CenA (Markowitz et 
al. 2007), MCG- 2-58-22 (Rivers et al. in press), Mrk 766 



(Turner et al. 2007), NGC 4258 (Reynolds et al. 2009), NGC 
4593 (Markowitz & Reeves 2009), NGC 5548 (Liu et al. 
2010) and NGC 7213 (Lobban et al. 2010). 



5 CONCLUSIONS 

Resulting from the work above on this small sample of five 
deep Suzaku observations of Seyfert 1 AGN, we conclude 
the following: 

(i) The warm absorber components (if present) in these 
AGN significantly influence the Fe K region, despite predom- 
inantly affecting the X-ray spectrum at lower energies. Fit- 
ting the full 0.6-100.0 keV broadband spectrum ensures the 
additional spectral curvature due to the warm absorber is 
accounted for prior to making any conclusions about the ac- 
cretion disc or central supermassive black hole. Indeed the 
broad 'red-wing' observed in AGN spectra appears much re- 
duced after a complete and consistent modelling of the soft 
X-ray spectrum. 

(ii) There is evidence for absorption features in the FeK 
region in 3/5 objects in this sample due to highly ionized 
zones of gas. Absorption due to FeXXV is found in MCG- 
6-30-15, NGC 3783 and NGC 4051 and FeXXVI absorption 
is found in MCG-6-30-15 and NGC 4051. The absorption 
found in NGC 4051 is blue-shifted (as noted by Lobban et 
al. 2011) and in MCG-6-30-15 (Young et al. 2005). 

(iii) Examining the FeK region after modelling the soft 
X-ray absorber and any highly ionized FeK absorption, we 
find the initially observed broad component in the data is 
significantly weaker. It is fitted well by kerrdisk line emis- 
sion from the inner regions of the accretion disc of mean 
strength EW = 99 ± 28 eV. Even in the case of MCG-6-30- 
15, the strength of the broad line is only EW = I6II44 eV, 
compared to a cumulative EW — 320^45 eV in Miniutti et 
al. (2007) in which the absorption zones are not fully con- 
sidered in a fit to the 3-45 keV spectrum. 

(iv) The use of a high column density {Nn ~ 10^^ cm"'^) 
partial covering component to account for the hard excess 
has only a marginal affect upon the derived accretion disc 
parameters. Subtle curvature is added to the spectrum at 
FeK energies due to the partial coverer, however, this is not 
enough to substantially alter the observed strength of any 
broad component which may be present. This is particularly 
notable in MCG-6-30-15, although this is not to say that 
some complex partial covering models cannot fully account 
for all the residuals in the FeK region. 

(v) Using the above conclusions, we estimate an average 
emissivity index of q — 2.8±0.2 in this small sample of AGN, 
indicating that emission is not particularly centrally concen- 
trated with the majority of emission from the accretion disc 
not occurring in regions approaching i?isco- Combining this 
with the results of Patrick et al. (2011), we find an average 
emissivity of g = 2.5 ± 0.2 over a sample of 10 objects. 

(vi) This analysis rules out a maximally rotating black 
hole with the kerrdisk model in Fairall 9, MCG-6-30-15 
and NGC 3783 at high confidence levels with increases in 
Ax^ of -1-6, -1-70 and +15 respectively. Fitting without a 
partial covering component similarly rules out a maximally 
rotating black hole where the fit is worse by Ax^ = +24 m 
MCG-6-30-15. 
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(vii) Given a full broadband model, including the re- 
quired absorption components, we estimate a BH spin of 
a = 0.49lo:?§ in MCG-6-30-15 at an emissivity oi q = 
2-71q^, accretion disc inclination of i° = 44^2 ''■nd a mod- 
erate strength of EW = IQltfi eV. 

(viii) Similarly we estimate an upper limit of a < 0.32 
and a < 0.94 for NGC 3783 and NGC 4051 respectively. 
The spin of the central BH in Fairall 9 is estimated to be 
a = 0.67lg'5^° whereas the spin cannot be constrained in 
NGC 3516. ' 

(ix) In NGC 3783 we find a marginal indication of ret- 
rograde spin, using the relline model yields a retrograde 
BH spin of a < —0.04 at the 90% confidence level with 
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Detailed tables of spectra parameters associated with fits 
performed in Section 3.2 and Section 3.4: 
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Table Al. Initial paraiiictrisatioii of the 0.6-100.0 kcV spectrum, modelled using powerlaw, warm absorption components, COMPTT, 
PEXRAV and Gaussians (all where appropriate). " POWERLAW normalization given in units (lO^'^ phkeV"^ cm~^ s^^). * Flux for cOMPTT 
quoted over the 0.6-lO.OkeV range in units 10~^^ergcm~^ s~^. Flux given in units (10~^ ph cm~^ s~^). Column density measured 
in units 10^^ cni^"^. Ionization parameter given in units ergcms"^. * Denotes a frozen parameter, f Emission line has been confined 
to the range 6.63-6.70 keV. 
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Table A2. Components for the baseline model (without broad or diskline emission) to the long observations with Suzaku XIS, HXD 
and BAT data from Swift. ° POWERLAW normalization given in miits (10^'^ phkeV^-*- cm^'^ s^'-). ^ Colmnn density for partial coverer 
measured in units lO'^'^cm^'^. Ionization parameter given in units ergcms^-'^. Flux for COMPTT quoted over the 0.6-10.0 keV range 
and REFLIONX over the 2.0-lOO.OkeV range in units lO^^^erg cm^"^ s^^ . '' REFLIONX normalisation given in units 10^^. ^ Flux given in 
units (10~® phcm~^ s~^). ^ Column density measured in units lO^^cm"^. * indicates a frozen parameter. The highly ionized zone for 
NGC 3783 has redshift fixed at z = 0.0097. 
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34.07l^;f* 

2 4+0.1 
^•^-0.1 

~ 50% 


06+0 04 

'^•^°-0.04 

q -1+0.5 
•^■1-0.5 

~21% 


- 


2 6+0-3 
~ 37% 


compTT 


kT (keV) 

T 

Flux'' 


<4.1 
-413 


< 8.8 

72+0-18 
"•'^-0.10 

-74 


- 


8 5+°-2 
°-°-0.1 

1 6+0-1 

9 24+5-31 

-461 


5 4+3-1 
^•^-1.0 

9+0-4 
u.y_Q 2 

49+0-08 
-333 


REFLIONX 


Norm" 
Flux'' 


< 11 

2 0+1-5 
^- "-'-0.35 


< 12 
1.0* 

1 02+0-09 
1 28+0-11 


< 12 

1.0* 

1 47+0-23 
9 no + 0.37 


< 14 

O.SlO;! 

5 oq+l-44 
0.0»_Q 26 

8 Qn+2-38 


< 21 

0.51^4 
0.86l«:g 
1 20+O-Sf' 

l-^'J-0.58 


FeXXV emission 


LineE (keV) 
EW (eV) 
Flux-'' 


fi fio+0.03 

23i« 
-40 


- 


6.6ltH? 

25+0-1'' 
-8 


- 


6 59+0-05 

16l^ 
0.49l°-i° 
-21 


FeXXVI emission 


LineE (keV) 
EW (eV) 
Flux-'' 

A 9 


7 01+0-05 
'•Ui_o.o5 

17+8 

-^■'-8 
'■'■^'■'-0.18 

-14 


6.97* 
lOtl 

-7 


- 


6 96+0-01 
44+11 

4 90+1 1^ 

-45 


- 


Warm absorption zone 1 


Nil 3 




r\ r»o + 0.()6 

0-22Io.03 

0.76+E!i» 
— u.uy 


1 nn+0.66 

l-90Io.2i 

— U.Ub 


^ /-O+0.03 

2.10+0-" 

— U.Ul 


r\ 1 + 2.09 

O-lS-o.ii 
3.01+0-3^ 

— U.o ( 


Warm absorption zone 2 


AThS 
log(0" 




0.47+°°8 

— U.Uo 




— U.IU 

0.69t°;lJ 


0.18+°-°^ 

— U.ii 

1 ao+0.07 


High 5 zone 


log«) " 

Vout (kms^-"^) 




3.99l?:l 
3200t«o 




4 Qn+0.80 
3 99+0-07 


q 55+18,95 

4.111!!:?;^ 

5600l«oo 


WABS 




0.0316* 


0.0392* 


0.178t°:°li 


0.0991* 


0.0115* 


XIS/BAT 


Const 


U.Dy_Q gg 


n Q^+o-06 


-1 09+0.11 

^•°^-0.08 


0.8410:03 


0.68l°:°| 






959.1/885 


2026.8/1823 


521.2/495 


1493.1/1379 


1091.5/1086 
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Table A3. List of soft X-ray emission and absorption lines included within the baseline model. " Flux given in units (10 ^ phcm ^ s 
negative flux and EW indicate an absorption line. 



LineE (keV) 


EW (eV) 


Flux" 




Identification 


0.66l°;« 




Fairall £ 
42 07+1*'-" 


1 

-18 


OVIII Lya 


0.681H; 
77+0.01 

'-0.01 

n QK+o.oi 
9 ,7+0.01 

„ -^+0.04 

2-77Io.o8 


411 

-Otl) 

-(2_2) 


MCG-6-30-15 
31.721^;^^ -31 

71.38tj3;o^ -892 
5-58tJ|^ -70 
-(2-65t°Ji) -196 
-(0.45l°:^e) -9 


OVIII Lya 
OVII RRC 

NelX resonance 
S XV ls-2p 
SXVI ls-2p 


o.solHi 

"•'^ -0.23 
1 44+0-02 

1 r;fj+0.02 


,7+10 

9+10 

5li 
8tl 


NGC 3516 
11.18l«;f, -7 
3.961^-35 

0.781°:^^ -11 
1.03l°j^ -21 


OVII RRC / FeXVII 3d-2p 
NelXHea 
MgXII Lya 
Mg XI He/9 


0.63l«;^J 
70+0 01 

'^^_0.01 

-, o/l+O.Ol 

1.47* 

Q 9Q + 0.03 

"J-^o_0.02 


10+2 

171^ 

0+1 
2-1 

i3i; 

-(131^2) 


NGC 3783 
197.10l^i:J» -30 

97.03l^^g^^ -146 

5.56i^;^2 

15.65ti;29 -359 
-(5.6010-^1) -163 
1.55l°;t5 -26 


VIII Lya 
OVII RRC 
JNeiX KKC / NeX ls-3p 
MgXII Lya 
SXV ls-2p 
Ar XVIII Lya 


9 99 + 0.03 
^•^^-0.03 




NGC 4051 
0.63«;« -16 


Si XII ls-3p 
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Table A4. Components for the dual reflector fit to the long observations with Suzaku XIS, HXD and BAT data from Swift. ° POWERLAW 
normalization given in units (10"'^ phkeV"^ cm"^ s~^). ^ Flux for COMPTT quoted over the 0.6-lO.OkeV range and REFLIONX over the 
2.0-lOO.OkeV range in units 10~^^ergcm~2 g— l c Ionization parameter given in units ergcms"'^. REFLIONX normalisation given in 
units 10~5. " Flux given in units (10 ^ phcm ^ s Column density measured in units lO^^ cm ^ . * Denotes a frozen parameter, f 

Emission line has been confined to the range 6.63-6.70 keV. Spin can only be constrained in MCG-6-30-15 and NGC 3783 at the 90% 
confidence level. The highly ionized zone for NGC 3783 has redshift fixed at z = 0.0097 and Vaut is fixed at the best fit value from 
the baseline model (TableEU for MCG-6-30-15 and NGC 4051. The COMPTT parameters for MCG-6-30-15 have been frozem at their 
best-fit values from the baseline model. 



Component 


Parameter 


Fairall 9 


MCG-6-30-15 


NGC 3516 


NGC 3783 


NGC 4051 


Powcrlaw 


r 


^■^'-'-0.01 


^•^■^-0.01 


1 76+0-04 
^^•'"-0.03 


1 oc+O.OS 
t-»5_0.01 


2 05+0 02 




Norm " 


n 07+0. 03 

'J-*''-0.03 


1 70+0.02 

^•'•^-0.02 


"■■^•^-0.01 


1 54+001 
^•^^-0.04 


03+0 05 




kT (keV) 


< 14.1 


3.9* 




> 9.5 


c: 7+8.0 
3- '-2.9 


compTT 


T 


U.O_Q 2 


0.8* 




1 g+0.1 






Flux*" 


'-'■■^"-0.02 


0.72* 




2 59+0-" 


49+0-04 






-449 






-2062 


-251 







< 58 


< 11 


< 20 


< 11 


< 30 


Unblurred REFLIONX 




2.0* 


1.0* 


1.0* 


1.0* 


1.0* 




Norm'' 


0.98tJ;g 


70+0-08 

"■'"-0.10 


72+0-21 
"■'■^-0.4,3 


4 50+'^-i2 


qg+0-10 




Flux'' 


2 42+4-28 
^■^^-1.34 


Q 7, +0.08 

'^•'^-o.io 




gQ9 + 0.19 


1 oo+"-i® 

'^■""-0.65 




Q 


7+2.5 
^- ' -0.2 


9 O+0.2 
^■•^-0.1 




9 0+0.4 






a 




o.6il«:i| 




< 0.31 




KERRCONV*REFLIONX 


i° 


33+5 


35t^ 




< 13 










< 11 


< 24 


52+^ 


< 18 






2.0* 


1.0* 


1.0* 


1.0* 


1.0* 




Norm'' 


26+°-^^ 
"■■^"-0.16 


2 42+0-28 
^•^^-0.27 


44+0-75 
'^•^^-0.18 




1 30+''l* 

1.00_g 72 




Flux'' 


n 79+2.06 
'J- '^-0.44 


9 90+0.26 
^•^°-0.25 


1 9Q+2.18 
J^-^°-0.52 


9 qi+0-96 


1 00+''" 

i.UU_g 7g 




LineE (keV) 




6.63+0-04 ^ 








Fe XXV emission 


EW (eV) 




g+5 
"-3 










Flux-^ 


43+0-30 
'-'■^■^-0.22 


25+0-20 
^^•^^-0.13 












-14 


-7 









LineE (keV) 7.021° 6.97* - 6.961^;! 

Fe XX VI emission W (eV) 20t^ 45l?| - 44lf 

Flux-^ 0.47^0.17 2.8ll^;°* - 3.40lg;' 

Ax^ -14 -10 - -49 



Warm absorption zone 1 



log(5) 



0.69t«:i? 



2.58 
2.09 



+0.39 
0.24 
0.05 
0.02 



4.61 
2.09 



+0.05 
0.02 
0.01 
0.01 



0.24 
2.77 



+0.31 
0.15 
+0.23 
0.44 



Warm absorption zone 2 



log(5) ' 



0.73 
1.77" 



+0.05 
0.06 
+0.04 
0.06 



0.15 



+0.24 



0.08 

< 0.86 



0.19 
1 



+0.11 

0.14 
+0.06 
-0.07 



^H^- - 11.0111^,0 - 1.59t«;i 13.2llJ«2 

\og{0- - 4.08l°:°9 - 3.87l°:°« 4.1llg:2i 

t)out(kms-l) - 3200* - - 5600* 



WABS 




0.0316* 


0.0392* 


n 1S9+0-013 
U.i»^_0.013 


0.0991* 


0.0115* 


XIS/BAT 


Const 


Q 7-, +0.06 
'J- '-'--0.05 


n SA+O-06 
U-«4_o,o5 


1-84_Q p7 


n. 04+0.02 


r. 77+0.05 
"■ ' '-0.06 






929.3/881 


2061.4/1823 


510.9/495 


1444.8/1378 


1111.2/1084 
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